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ABSTRACT 



We propose to use gravitational wakes as a direct observational probe of the 
collisional nature of the dark matter. We calculate analytically the structure of a 
cn i wake generated by the motion of a galaxy in the core of an X-ray cluster for dark 

matter in the highly-coUisional and coUisionless limits. We show that the differ- 
^ i ence between these limits can be recovered from detailed X-ray or weak lensing 

v^ i observations. We also discuss the sizes of sub-halos in these limits. Prelimi- 

O ! fluid-like dark matter but does not exclude scenarios in which the dark matter is 

weakly collisional. 

O ' Subject headings: dark matter - X-rays: galaxies: clusters - galaxies: clusters: 

-t— » ' general - hydrodynamics 

> 

/\ ' 1. Introduction 



Although the dark matter is known to dominate the mass budgets of bound systems 
such as galaxies or galaxy clusters, its nature remains mysterious. On cosmological scales 
(^ 1 Mpc), the standard "cold dark matter" (CDM) model, in which the dark matter 
is non-baryonic, non-relativistic, and coUisionless, has had considerable success matching 
observations. However, two sets of problems have emerged in applying the same model to 
smaller scales ( ^pergel fc Steinhardt 2000| ; |Hogan fc Dalcanton 2000] and references therein) . 



The first is known as the "substructure problem." The CDM model predicts significantly 
more substructure (in the form of relic dark matter halos) on small scales than is observed. 
For example, models of the formation of the Local Group of galaxies predict an order of 
magnitude more dwarf halos than the observed number of dwarf galaxies ([Moore et al. 



|1999| ). The second problem is known as the "cusp problem": with no interactions to prevent 
collapse, CDM halos should form nearly singular cores (e.g., [Navarro, Frenk, fc White 1997| ). 
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In contrast, observations indicate that galaxies have flat dark matter profiles near their 
centers ( [I'lores fc Primack 1994| ; Borriello & Salucci 2001). 

These problems led Spergel & Steinhardt (2000) to reintroduce self-interacting dark 



matter (SIDM), originally proposed by Carlson, Machacek, & Hall (1992; see also |de Laix 



Scherrer, fc Schaefer 1995| ). In this scheme, the dark matter has a significant cross-section 



for self-interaction (and possibly a significant cross-section for interaction with baryons as 
well; see [Kusenko fc Steinhardt 2001| ). In such a picture, cusp formation is halted by heat 



conduction (which heats up the dark matter at the halo centers) and the formation of shocks. 
The substructure problem is avoided because dark matter particle interactions evaporate 
small halos on relatively short timescales. 

Existing limits from laboratory experiments restrict the range of cross-sections available 
to SIDM that interacts directly with baryons ( pVIcGuire fc Steinhardt 200 1|) but place no 



restrictions on the cross-section for dark matter-dark matter interactions. To constrain this 
cross-section, we must turn to astronomical observations. In this paper, we propose such a 
test, based on observations of the motion of individual galaxies inside galaxy clusters. 

As a galaxy moves through a cluster, the gravitational potential of the galaxy induces 
a wake in the surrounding medium. The back-reaction of this wake on the galaxy is called 
dynamical friction and is known to be important in the dynamics of cluster galaxies (as well 
as satellite galaxies and globular clusters on smaller scales; pinney fc Tremaine 1987| , [Frenk 



et al. 1996]) . Wakes are induced in both the dark and gaseous components of the intracluster 



medium (ICM). The structure of the dark matter wake provides a signature of the collisional 
nature of the dark matter that is potentially observable through its gravitational lensing 
effect. The total gravitational potential also imprints a wake in the gas distribution; a precise 
measurement of the morphology of this wake yields the size of the perturber's dark matter 
halo. This in turn constrains the collisional nature of the dark matter: the sizes of CDM 
halos are set by tidal stripping in the cluster potential ( [Merritt 1988|) , while ram pressure 



stripping truncates SIDM halos at considerably smaller radii. The exquisite resolution of 
the Chandra satellite allows very precise measurements of the extent of surface brightness 
variations in the ICM gas and could potentially constrain the collisional properties of the 
dark matter. 

In order to approach this problem analytically, we will focus on two extreme scenarios 
for the nature of the dark matter. The first is the CDM model, in which the dark matter 
is collisionless. The second model assumes that the dark matter is strongly self-interacting, 
so that it may be treated as a perfect fluid. We will refer to this scenario as the Fluid 
Dark Matter (FDM) model. The latter approximation has been used before to isolate the 
implications of SIDM models ([Yoshida et al. 2000| ; [Moore et al. 2000| ). Our calculations will 
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thus bracket the parameter space available for dark matter self-interaction cross-sections. 
Our experimental test is easiest to interpret in the limit where the perturber mass is much 
smaller than that of the host cluster, and so we use linear theory in our FDM calculations. 

We first outline our formalism for calculating the wake structure induced by a perturber 
moving through coUisionless and coUisional media in §2. We then describe the effects of the 
ICM dark matter on the perturber halo in §3. Finally, we present our results and compare 
them to simulations in §4, and we discuss their implications in §5. 



2. Gravitational Wakes Due to Dynamical Friction 

First, we develop the mathematical formalism necessary to calculate the density distur- 
bance in a cluster due to the gravitational effect of a moving galaxy or group of galaxies. In 
§2.1 and §2.2, we calculate the wake in coUisionless and coUisional media, respectively. In 
a coUisional medium, an additional wake will be created through deflection of the ambient 
medium around the surface of the perturber. We discuss this "hydrodynamic wake" in §2.3. 
We assume throughout that the perturber (e.g., a galaxy) with mass Mg moves through a 
stationary, uniform background medium (e.g., a cluster core) of density pc with a constant 
velocity Vg = M.ac directed along the 2;-axis, where M. is the Mach number of the flow 
and (Jc is the velocity dispersion of the background medium (assumed to be equal to the 
sound speed for simplicity). We include a discussion of when the requirement of a uniform 
background density may be relaxed in §4. 



2.1. CoUisionless Medium 

First we consider the structure of the wake induced in a coUisionless medium, following 
the method first introduced by Mulder (1983; see also [Weinberg 1986|) . Without collisions. 



the particle distribution function /(x, v) must be a function only of the integrals of motion 
( [Binney fc 'Iremaine IQS'^ ), which can be calculated for particles deflected by the perturber 
using planar orbit theory (Mulder 1983). The precise form of the distribution function 
is then determined by the boundary conditions: in the frame of the perturber, / must 
approach a Maxwellian distribution with mean velocity — Vg and velocity dispersion ac at 
infinity. Finally, the density is p = / fd^^v. Note that, given the assumption of a uniform 
background medium, this treatment is exact, and it works equally well for a point mass or 
an extended body. 

To linear order (that is, neglecting the small fraction of background particles that are 
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strongly deflected by the perturber), the wake due to a point mass can be treated analytically 
( [Mulder 1983| ). In the frame of the perturber, the density enhancement A = (p — l)/p is 



AcL,pt(x, t) ^ — - — exp — (1 - cos I 



cr2 r 



fMcosd 



where (r, 6) are polar coordinates centered on the perturber's current location. (Note that 
Mulder's derivation contains a sign error, though his final expression is correct. We find 
e ^ —6 in Mulder's equations [AIV.3] and [AIV.4].) This approximation can be derived in 
a completely independent manner using the theory of linear response ( [(Jolpi fc Pallavicini| 



19981) . Note that the density structure is independent of the propertiers of the perturber 
and the background medium if we rescale r = t/Rq, where Rq = GMg/a"^; this is true for an 
extended perturber as well. 

Because of the lack of shocks, wakes in a coUisionless medium are smooth regardless 
of the velocity of the perturber, with the peak density enhancement directly behind the 
perturber. Increasing the velocity focuses relatively more of the wake behind the perturber. 



forming a center- filled cone-like structure at high velocities (see [Mulder 1983| , Fig. 2). This 



is to be contrasted with the well-defined, edge-filled cones we find below in collisional media. 



2.2. Collisional (Fluid-like) Medium 

We calculate the wake in a highly collisional, or fluid, medium using linear perturbation 
theory. The passage of the object perturbs the background medium, yielding a density field 
p = pc(l + Ap); we assume that Ap ^ 1. The first-order continuity equation is then (see 
also Pstriker"T999D 

1 HAr. 

+ V ■ vi = 0, (2) 



1 dAp 



ac dt 
where vi is the perturbation to the velocity field, and the Euler equation is 

^^'''+VAF = -lv$ext, (3) 



(Tc dt a^ 

where $ext is the gravitational potential of the perturbing body. Combining these equations 

with Poisson's equation for the potential, eliminating vi, and Fourier transforming over both 

space and time yields 

/,2 cu2\ . 47rG. 

Here Pg = Pg{R,z — Vgt) is the specified density field of the perturber and a hat denotes 
a Fourier transformed quantity. This may be simplified by noting that Pg = (27r)^/^(5(u; — 
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kzVg)po, where po is the density distribution of the perturber held stationary at the origin. 
Taking the inverse Fourier transform and evaluating the Dirac delta function, we find 

Ap(x.i) = ^/rf^k ,, ^°„., e--. (5) 



where s = x — f^tz. The integral on the right may be inverted with the convolution theorem. 
However, the inverse Fourier transform of /(k) = l/{k'^ — Ai'^k'^) is not trivial ii Ai > 1, 
because in that case the integrand has poles at k^ = ±{kl + kyY^'^/P, where (3"^ = M.^ — 1. 
We therefore evaluate the integral by transforming to the complex plane and choosing a 
contour slightly above the real axis. If Sz > 0, the point of interest is ahead of the perturber 
in real space, and causality dictates that the medium remains unperturbed. We close the 
contour at +ioo in this case. Conversely, if s^ < 0, we close the contour at —ioo. Then, 
using the convolution theorem, 

Ap(x,t) = ^ fdh'- ^-^P^ -^, (6) 

^c J (s, - z'Y - i3\R^ + i?'2 - 2RR' cos Q'^^ 

where 



2 A^ > 1, (s^ - /) + /3i/(i?2 + i?'2 - IRE! cos d') > 0, 
e'=<| 1 -M<1, (7) 

otherwise. 

In the case of a point mass, po = Mg6{r'), this simplifies to the result found by Ostriker 
(1999): 



^F,pt 

where 



AF,pt(x,t) 



^iz-Vgty-f3^R^'' 



2 M> l,{z-Vgt) < -f3R, 
C={ I M<1, (9) 

otherwise. 

As in the collisionless case, the wake structure is independent of the properties of the per- 
turber and the background medium if we rescale s = s/Rq. Note that these formulae are not 
valid in the case TW = 1. As shown by Ostriker (1999), the dynamical friction force due to 
the wake increases rapidly as the velocity approaches the sound speed, so such a perturber 
quickly slows to subsonic speeds. 

In the collisional case, the structure of the wake depends strongly on whether the motion 
is subsonic or supersonic. In the subsonic case, the wake has front/back symmetry relative 
to the perturber and varies smoothly. However, in the supersonic case, the wake forms a 
Mach cone trailing the perturber with opening angle cot~^ (3. The overdensity at the edges 
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of the Mach cone is infinite for a point mass. An extended perturber smootlis tlie singularity, 
spreading tlie mass contained in tlie edge over a region of width ~ 2Rg. 

For a point mass wake, the column density 5S along any sightline intersecting the cone 
can be computed analytically and is independent of the location in the wake, 

where i is the angle between the direction of motion of the perturber and the line of sight. 
In the case of an extended perturber, the above formula is very close to the true column 
density except near the edges of the cone, where the column density drops to zero over a 
width ~ 2Rg. 

The singular edge to the cone in the point mass case, and the concentration of mass 
in this region for an extended mass may raise some concern about the linear perturbation 
analysis used here. Simulations have shown that, even for a point mass, the resulting dy- 
namical friction force is nevertheless quite accurate ( ^anchez-Salcedo fc Brandenburg 1999| ) , 
so we expect the linear treatment to work adequately. We do find that Ap < 1 throughout 
most of the volume of interest (far from the perturber), so linear theory should be accurate 
in these regions. We include a more in-depth discussion of existing simulations in §4.3. 



2.3. Hydro dynamic Wake 

For a collisional medium, we must add the hydrodynamic disturbance introduced by the 
passage of the perturber to the gravitational wake. If the bound halo is "stiff," it deflects 
ICM particles and sends pressure perturbations into the ambient medium. The hard surface 
approximation is expected to be adequate for the ICM gas, because the Larmor radius for 
typical cluster magnetic fields is very small (^ 1 pc), although particles may still penetrate 
the perturber if the field lines intersect the body. The width of the edge of the dark matter 
halo is expected to be of order the dark matter particle mean free path. 

The density disturbance due to the surface effect may be calculated analytically far 
from the perturber if the perturbing body is streamlined ( [Landau fc Lifshitz 1959( , §115). 
We switch to the rest frame of the perturber and use coordinates s with the origin at the 
center of the body. In this frame, the gas velocity is vq = — foSz + vi and the density is 
p = Pc(l + Ahd), where Ahd is the density perturbation due to this hydrodynamic effect. 
Far from the perturbing galaxy, the gas is only slightly disturbed and linear theory applies. 
If we assume that vi is related to a velocity potential (p via vi = V0, the Euler equation 
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inay be written as 

where we have assumed that and Ahd vanish at infinity. Using the continuity equation to 
ehminate the density enhancement, we find 

dsl ds^ dsl 



The solution is ( [Landau fc Lifshitz 1959|) 



Vq fi^ S'{u)du 



2vr Jc ^(s^ _ uY - /?2i?2 

where C = max(s2 + PR, —i), 2i is the length of the perturbing body in the direction of 
motion, S{z) is the cross-section of the perturbing body, and a prime represents differentia- 
tion with respect to u. The length i of the body is not well-determined in our analysis; the 
linear theory is strictly only valid for bodies with very small opening angles. We let i = xRg, 
where Rg is the transverse radius of the object (see §3). Note, however, that the assumption 
of a streamlined perturber is not likely to be valid in realistic circumstances. 

In this limit, the hydrodynamic wake forms a cone with a structure similar to the gravi- 
tational wake. Both are "edge-filled" cones with opening angle cot~^ f3. In the approximation 
used above, weak discontinuities appear at the Mach cones originating from the leading and 
trailing edges of the body; in reality, far from the perturber, these are weak shocks ( [Landau 



fc Lifshitz 1959[ ; §114). The ICM gas is compressed at the leading shock and then becomes 



rarefied as it defiects around the body. The trailing shock discontinuously increases the 
gas density. As equation ( [TB| ) shows, the density distribution of the wake depends on the 
cross-section of the body, which is not well-determined by our treatment (it depends on the 
detailed initial mass distribution, gas replenishment in the galaxy, and the nature of the 
ram-pressure stripping). Unfortunately, the dependence on the cross-section is substantial, 
particularly between the two discontinuities. In addition, real halos are not streamlined 
( [Markevitch et al. 20001 ; [Vikhlinin, Markevitch, fc Murray 2001[) , but rather have blunt sur- 



faces. In this case the leading shock is "detached" from the front of the halo. The bow shock 
is normal to the direction of motion on the collision axis and asymptotically approaches the 
Mach cone at large distances ( [Schreier 1982[ ). Within this region, the hydrodynamic motions 



are quite complex, but some general conclusions can be drawn. 

Behind the trailing shock, the density enhancement due to the perturbing body falls off 
rapidly. This is best illustrated by evaluating equation (p!3[) for i? <^ — s^. We find 

AuB^VxM'l^] , (14) 



where 77 ~ 1 depends on the cross-section of the perturber. We find that Ahd ^ 1 throughout 
most of the cone interior. However, the compressions and rarefactions are strongest in the 
region bounded by the two shocks discussed above. We find that in this region 



Ahb-^^ m^ P[ ^^ -) --R.<s. + PR<xR,. (15) 

Here p{y) is a function that depends on the precise form of the perturber's cross-section; in 
general, it changes sign as R varies. Because Sz + /3R is bounded, in this region 5Shd oc 
Ahd oc |sz|~^/^. We find for some generic cross-sections that |Ahd| > 1 within the region 
of interest, indicating that hnear theory is not vahd. Therefore we cannot meaningfully 
compare the calculated amplitudes of the variations in this region to the gravitational wake 
without detailed simulations of the interaction process. However, note that the integrated 
mass perturbation along any line of sight passing near the collision axis 5Shd vanishes. As 
the body moves, it produces finite outgoing cylindrical sound waves. The total displaced 
mass through any line of sight must therefore vanish ( [Landau fc Lifshitz 1959| , §95). 

There is yet another effect that can create a wake in a coUisional medium. Stevens, 
Acreman, & Ponman (1999) simulated the effects of ram pressure stripping on the interstellar 
medium of galaxies moving through the ICM. Most of the gas is stripped into a plume trailing 
the galaxy, but they also found "bow shock" wakes with similar structure to those produced 
by point masses. The simulation box size was comparable to the size of the dark matter 
distribution of the galaxy, and so it is not possible to infer directly the effects of the stripped 
gas on large scales (^ -Rgas) which are of interest here. The finite mass of gas in the 
galaxy (particularly after it has traveled some distance through the cluster) indicates that 
the stripped gas should not have an appreciable effect on the large-scale wake structure (see 
§3.2 below). 

In summary, we expect that hydrodynamic interactions around the perturber will induce 
a wake of similar structure to the gravitational wake in a coUisional medium. Although the 
overdensity in the interior of the hydrodynamic wake falls off quickly with distance from 
the perturber, the mass between the shocks produced by the leading and trailing edges 
of the pertuber decreases much more slowly and, in fact, linear theory breaks down for a 
large region around the perturber. Therefore we cannot put meaningful constraints on the 
observability of the hydrodynamic wake. 



Halo Truncation in Galaxy Clusters 
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3.1. Dark Matter Stripping Processes 

Next we describe how to apply the above formahsm to galaxies and groups moving 
through X-ray clusters. We assume that the perturber moves through a constant density 
cluster core, surrounded by a softened isothermal envelope 

""^"^ = 1 2nJ[l + {r/r^r] r > Z ^^^^ 

where Uc is the total number density of ions and electrons in the core and Tc is the core 
radius. For simplicity, we assume the dark matter distribution in the perturber to be a 
singular isothermal sphere with a specified velocity dispersion cr^, so that Pg = 0"|/(27rG'r^). 
The mass contained within a radius Rg is then Mg = 2a^Rg/G. 

The radius of the halo is determined by the nature of the dark matter. If the dark 
matter is coUisionless, only gravitational variations within the cluster can affect the halo. 
Thus tidal stripping determines the final size, and the truncation radius may be found by 
solving Pg{< Rg) = Pc, where Pg{< Rg) is the mean density inside Rg ( [Babul fc 'laylor 2000] ), 
as also implied by numerical simulations ( [Merritt 1988| ) . If the dark matter is strongly self- 
interacting, the halo is subject to ram pressure stripping as well. The truncation radius is 
then determined approximately by balancing the ram pressure and the external pressure, 
Pcv'^ = pg{Rg)a'^. This simple prescription ignores two competing effects. First, the super- 
sonic motion of the perturber creates a bow shock leading it (see §2.3), which shields the 
halo from the ICM ram pressure. However, simulations show that for mildly supersonic 
motion the bow shock provides relatively little protection ( [Stevens et al. 1999| ; [Quilis et al 



2000 ). Instead, recent high- resolution simulations of galaxy motion through the ICM show 



that turbulent and viscous stripping are as effective as ram pressure in truncating gas disks 
( [Cjuilis et al. 20ll(]| ). More realistic prescriptions will therefore only increase the contrast 
between the FDM and CDM model predictions. Solving for Rg in the two cases yields 

-3 „_-3 f n «\ 1/2 



5x 10"^ cm"^ fg 0.6 V ( ^a 



where 



^'"^''H ^. ^- l300k^)"™''P'=- (^^) 



_J (a/37W) ^{ag/ac), coUisional medium, , , 

I 1, coUisionless medium, 



fg is the gas-to-dark matter mass ratio in the cluster, p is the mean atomic weight per 
particle, and the subscript "DM" refers to the dark matter medium. The length scale Rq is 



n^ 0.2 p I \ 0.3 / V300 km s 
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Note that the scaled size of the perturber is 

Studies of galaxies in nearby clusters, particularly Virgo and Coma, have shown con 
vincing evidence for ram pressure stripping of the interstellar atomic gas ( |Warmels 1988 




Uayatte et al. 1994| ). Magri et al. (1988) showed that the HI flux from spiral galaxies in- 



creases with distance from the cluster center and decreases with the X-ray luminosity of the 
host cluster, as one would expect for ram pressure stripping. Most recently, Bravo-Alfaro 
et al. (2000) examined 19 bright spiral galaxies in the Coma cluster and found a clear cor- 
relation between depth in the cluster and HI deficiency relative to isolated spirals. Rubin, 
Waterman, & Kenney (1999) reach a similar conclusion for the ionized gas disks of spirals 
in the Virgo cluster. 

It is much less clear whether dark matter may be stripped in these environments. Despite 
early evidence showing that cluster galaxies tend to have falling rotation curves ( |Whitmore,| 
[b'orbes, fc Rubin 1988| ), more recent studies have not found any systematic difference in the 
shape of the curves with location in the cluster ( pale et al. 200l| , and references therein). 
However, these observations do not rule out the FDM model, because the observed rotation 
curves are based on Ha emission. This gas is subject to ram pressure and is therefore 
at least as compact as the dark matter. A better test is to examine the stellar rotation 
curves, which at these distances can only be probed with bright tracers such as planetary 
nebulae ( |Ciardullo, Jacoby, fc Dejonghe 1993| ; Arnaboldi et al. 1996| ; [Arnaboldi et al. 1998| ) . 



To date, such studies have been able to examine the inner ~ 20 kpc of relatively massive 
ellipticals, about twice the extent of the emission line rotation curves. The observations show 
no evidence for deviations from the normal fiat rotation curves; however, their limited extent 
would not exclude a scenario that is only somewhat less efficient than our FDM model. 

Another limit comes from the fact that as the halo loses its dark matter, luminous 
matter in the galaxy becomes unbound. As a simple model, we consider a spherical galaxy 
with both the dark matter and stars distributed as singular isothermal spheres. In this case 
the distribution of stellar velocities is Maxwellian with standard deviation ag ( [Binney "S 



i'remaine 1987] ). We suppose that the initial radius of stellar matter is Ri. After stripping. 



the dark matter radius is reduced to Rg] if Rg < Ri, we assume that stars outside of Rg can 
escape from the galaxy if their total energy is positive. The fraction of stars that escape the 
galaxy is then 

/esc = ^ / " erfc I J- I dr. (21) 
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Using the observed relation between the velocity dispersion of a halo and its effective radius 
i?e ( [Binney fc IVlerrifield 1998| , p. 206 ), we find that in the FDM model 



^^IIM /^ 5xl0-^cm-^ /, 0.6\-^/^ / a, \-°-^ / ^c 

Rg^ ' V ^c 0.2 fx J V200kms-V UOOO km s" 

where we have set Re = 0.36i?j for a singular isothermal sphere. The well-known "funda- 
mental plane" relates the luminosity (assumed proportional to the number of stars), the 
effective radius, and the velocity dispersion of ellipticals. The scatter in the relation is only 
~ 0.07 dex ( Mobasher et al. 1999 ). Because the mass loss mechanism described above 



moves galaxies away from the fundamental plane, the scatter implies that Ri/Rg <^ 1.2 for 
galaxies in clusters, or that A4 ^ 1.1 in typical rich clusters. However, note that we have 
neglected the contribution of the (non-stripped) stars to the mass of the galaxy, which may 
be substantial in the inner regions. This limit is therefore currently no stronger than that 
provided by rotation curves of galaxies. 

On the other hand, observations of dwarf galaxies in the Local Group indicate that dark 
matter stripping may be much more effective than the CDM model predicts. At least three 
dwarf spheroidals (Carina, Draco, and Ursa Minor) are observed to have stellar tails outside 
of their optical cutoff radii ( |lrwin fc Hatzidimitriou 1995| ) . Significant amounts of dark matter 



outside of the optical radii would preclude escape of these stars, indicating that Rg ~ 300 pc 
for these systems ( [Moore 1996| ; [Burkert 1997| ). Our model predicts -Rg,cDM ~ 5 kpc and 
Rg,FDM ~ 200 pc for a typical dwarf galaxy with ag ~ 10 km s~^, favoring a model that is 
"nearly" FDM. One caveat is that, if the orbits of the dwarf galaxies are highly eccentric, 
tidal stripping may be much more effective than is inferred from their present locations. 

Heat conduction from the hot ICM to the bound halo dark matter will also tend to 
evaporate the halos. However, the thermal conductivity of a medium is k oc A, where A is 
the mean free path of the particles. In the extreme FDM limit, A = and no heat conduction 
occurs; Gnedin fc Ostriker (2000) find that heat conduction may be neglected over a Hubble 
time if a-Du/mDM ^ 10"^ cm^ g~^, where ctdm and moM are the self-interaction cross-section 
and mass of the dark matter particles, respectively. If A ^ 0, heat conduction and ram 
pressure will work together to truncate halos, making that process more efficient than we 
would otherwise expect. In the opposite regime of very weakly interacting dark matter 
particles, scattering interactions between halo and cluster particles can also evaporate the 
halo. However, in the CDM limit, this process can be neglected. Gnedin fc Ostriker (2000) 
find that the evaporation time is greater than a Hubble time for (Tdm/'^dm ^ 0.3 cm^ g~^. 
Note that evaporation operates independently of the motion of the galaxy through the ICM, 
and it will therefore truncate all subhalos, even those of central cD galaxies. 
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3.2. Indirect Effects of Dynamical Friction in Clusters 

A common difficulty in all treatments of dynamical friction in a uniform infinite medium 
is that the mass contained in the wake diverges with distance from the perturber. This 
divergence is normally avoided by assuming that the perturber is embedded in a medium 
of finite extent ([Binney fc Tremaine 198^ ). The divergence implies that there must exist 



^^"°^ 1, 0.3 ) 


'giM)-'/'. 


FDM, 


W2(l-^2) 


1- 


-eif(^Mfi/V2) 


rf/i, CDM. 



a distance r„ake inside of which the wake contains mass Mg. Beyond r^ake, we expect the 
second-order effect of the dark matter wake structure on the observed medium to dominate 
the effect induced by the perturber itself. Using the analytic approximations for the wake 
due to a point mass, we find 

fwake ^ ISa^J, ( ^^ 1 giM)-'/'. (23) 

where 

^^ ' \ I\ exp [-^M\l - fi^)] [l - erf (A^/x/v^)] rf/i, CDM. ^ ^ 

In the cases we study, we find that f„ake is comparable to or somewhat greater than the box 
size of our calculations (and thus larger than cluster cores; see §4). Therefore, we neglect the 
effect of the dark matter wake on itself and on the gas wake. We have confirmed numerically 
that the dark matter wake has only a small effect on its own structure within f„ake- We can 
also examine the extent of the stripped gas wake. If we let Mgas = kMg and assume that 
the stripped gas forms a wake similar to the one we have calculated, the stripped gas will 
be contained within a distance fgas ~ k^^'^f^ake -C r^ake- Thus the stripped gas should be 
confined to the central region around the perturber. 

Note that the dynamical friction force Fdf oc M^ in both the CDM and FDM models, 
with the coefficient /df varying by a factor of a few between the two ( Ostriker 1999 ). Because 



FDM halos are less massive than CDM halos by a factor of q;fdm, we expect that FDM halos 
will lose energy to dynamical friction at a slower rate than CDM halos and will therefore 
sink to the cluster core at a slower rate. However, FDM halos experience a ram pressure 
force Fram as well. In our model, -FoFZ-^ram ~ 4JDF/(A^^-Rg) -C 1, so the FDM sink time is 
determined primarily by ram pressure. The different causes of sinking in the two regimes 
can lead to observable effects in the galaxy distribution ( [Moore et al. 20001) . In particular, if 



ram pressure is effective, the time for a galaxy originally situated a radius Vf, from the cluster 
center to sink to the center is tram ~ 4=fb/<^c- This is generally much smaller than the sink 
time for CDM halos; we find 

^DF.CDM 4.4A^3 ^ g^/g \^ 



f^U 



t 



ram, FDM -'DP 
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Like Moore et al. (2000), we find that the sink time in the FDM regime is independent of 
the properties of the galaxy. In contrast, the sink time in the CDM regime has a fairly steep 
dependence on (Xc/cg, implying that in a given cluster we should observe the massive galaxies 
to be more dynamically relaxed than smaller galaxies. Such mass (or luminosity) segregation 
has been observed, although its origins are controversial because morphology segregation in 



clusters can mask intrinsic mass segregation ([Kashikawa et al. 1998| and references therein). 



Kashikawa et al. (1998) have found that brighter galaxies are more centrally concentrated 
than faint galaxies in Coma, though the degree of concentration depends on morphology. 
Furthermore, Drinkwater, Gregg, & Colless (2001) found that the velocity dispersion of the 
giant galaxies in the Fornax group is a factor of ^J2 smaller than the velocity dispersion of 
the dwarf galaxies, indicating that the former are in a more relaxed state. This may be due 
to the more rapid effects of dynamical friction on larger galaxies. A similar effect has been 
observed for galaxies in rich clusters ( [Adami, Biviano, fc Mazure 19981) . 



A further observable effect that we have not included in our calculations is the possibility 
of a separation between the dark and gaseous components of the perturber. In the FDM 
case, the two components are subject to the same decelarating forces (ram pressure and 
dynamical friction), but a CDM halo feels only dynamical friction. We therefore expect 
CDM halos to lead gaseous halos by an amount that depends on the time history of the 
interaction. Approximate equilibrium can be reached if the gravitational attraction between 
the dark matter halo and the gas halo balances the ram pressure experienced by the gas 
halo (note, however, that such a separation may be convectively unstable; [Ricker fc Saraziii 



20011) . We can estimate this distance if we ignore other forces (including gradients in the 



cluster potential); we then find that the equilibrium separation is 

/ 7 5 X 10-^ cm-^ O.ey/' / /, \ (oJaA'l^ ( a, \ 
^- = 'H:^ ^c ^) Mi^ Uokms-0 'P^' ^2'^ 

where 7 ~ 1 depends on the structure of the gas and dark matter halos. Such a separation 
has been observed between the galaxies in the NGC 4839 group and its intragroup gas as the 
group falls into the Coma cluster, and it can be attributed to the much smaller ram pressure 
experienced by the compact galaxies ( [Neumann et al. 2001[ ). The cluster parameters of 



Neumann et al. and the observed values of Ogjoc = 0.3 and A^ ~ 1.7 ( [Colless fc Dunn 1996[ ) 
suggest an equilibrium separation of D^q ~ 1 Mpc, much larger than the observed separation 
of ~ 200h^^ kpc, where the Hubble constant Hq = lOOh km s~^ Mpc^^. This is consistent 
with other suggestions that the merger is recent, with the NGC 4839 group on its first infall 
into Coma. Nevertheless, combining the X-ray map with a weak lensing map in this or a 
similar system could constrain the coUisional nature of the dark matter simply by locating 
the mass peak. 
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4. Results 

Two techniques may be employed to observe the effects of the truncated dark matter 
halos. First, gravitational lensing can constrain the projected surface mass density of the 
perturber/cluster system (§4.1). Alternatively, the wake in the ICM gas can be observed 
through its X-ray bremsstrahlung emission (§4.2). The former method has the advantage 
of directly testing the dark matter profile, but the X-ray signal is stronger and easier to 
observe. 

In presenting our results, we must ensure that our formalism (in particular the assump- 
tion of a constant density background medium) remains valid. If the background medium 
has Po = Po(-R)) where R is radial distance in the cluster, we have 
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in place of the wave equation with a single source term that led to equation (H). We expect 
that (dAp/dR) ~ a~^{d^^^t/dR) ~ I/Rq while po^dpo/dR) ~ l/r^. Therefore, in order 
to neglect the bracketed term on the right-hand side, we require Tc/Rq = fc 3> 1. A proper 
treatment in regions with rapidly varying density (as well as an examination of the effects 
of shocks and other nonlinearities) requires numerical simulations. For these reasons, we 
discuss our results in light of existing simulations in §4.3. 



4.1. Projected Surface Mass Density 



The top two panels of Figure 1 show the surface density of a a^ = 500 km s^^ group 
merging with the core of a cXc = 1000 km s~^ cluster with ric = 5.4 x 10~^ cm~^ and Vc = 
0.4 Mpc. The group has A^ = 2. The top left panel shows the surface density for the FDM 
model, and the top right panel shows the corresponding result for the CDM model. The 
horizontal dashed lines in each panel mark where the Mach cone with its apex at the center 
of the perturber intersects the edge of the cluster core; the apparent suppression of the signal 
below this line occurs only because we do not allow the wake to propagate out of the cluster 
core. 

The most straightforward way to detect the FDM Mach cone is to look for the abrupt 
drop in projected mass at its edge. The strength of this signal can be estimated by comparing 
the (constant) surface density of the cone, given by equation (||), to the surface density of 
the cluster in which it is embedded, S. In the limit in which the impact parameter of the 
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Fig. 1. — Results for a a^ = 500 km s~^ perturber moving with A^ = 2 through a background 
cluster with ac = 1000 km s~^, ric = 5.4 x 10~^ cm~^ and Vc = 0.4 Mpc. Top left: Projected 
mass profile, FDM model. Top right: Projected mass profile, CDM model. Bottom left: 
Surface brightness, FDM model. Bottom right: Surface brightness, CDM model. Contours 
are set at 1% of the peak surface density of the background cluster in the top panels and 
at 5% of the peak surface brightness of the background cluster in the bottom panels. The 
merger occurs in the plane of the sky with the perturber moving along the z-axis. Axes 
(in cylindrical coordinates) are in units of Rq appropriate to each panel (see Eq. ||19|| ). 
Horizontal dashed lines mark the intersection of the Mach cone with the cluster core. 
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line of sight is much less than r^, the excess is 
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The amphtude depends strongly on the ratio of the velocity dispersions, so the best can- 
didates are galaxies moving through groups or groups of galaxies moving through clusters. 
The small surface density of groups makes lensing observations difficult, so we focus on the 
latter possibility in this section. 

As discussed in §2.3, coUisional media also give rise to wakes from hydrodynamic and 
stripping processes. The effects of the hydrodynamic wake well behind the perturber should 
be small, because (5Shd vanishes to linear order. However, the stripped dark matter could 
enhance the signal immediately behind the galaxy and the leading bow shock may be quite 
strong as well; we therefore regard our calculation as a lower limit to the true enhancement 
in the FDM case. 

Another observable difference between the two pictures is the amplitude of the mass 
perturbation induced by gravity. Tidal stripping in the CDM model leaves halos much more 
extended than FDM halos. They are therefore more massive and have correspondingly more 
gravitational influence on the ambient medium; equation (EBI) shows that (5S/S oc odm- In 



fact, the enhancement in Figure 1 is somewhat smaller than this because the CDM wake 
does not have the same structure as the FDM wake. 

A third observable effect is the size of the dark matter halo itself, which is dramatically 
different in the two regimes because of the different stripping mechanisms. By correlating 
the mass maps (from lensing) with X-ray maps, we can compare the relative sizes of the 
gas halo (which we know to be stripped by ram pressure) and the dark matter halo. One 
method is that of Natarajan & Kneib (1997), who analyzed the variation in the shear of 
background galaxies in concentric annuli around cluster galaxies in order to constrain the 
sizes of their dark matter halos. So far, the method has yielded only an upper limit on the 
sizes of the halos that is comparable to their tidal radius ( Natarajan et al. 1998 ). Such a 



comparison should be more straightforward in the case of a group of galaxies merging with 
a larger cluster. Note that this effect does not necessarily rely on the extreme FDM model, 
because a combination of ram pressure, heat conduction, and evaporation due to scattering 
will occur across the entire possible range of parameters for SIDM. 
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The recent Chandra images of sharply defined gas halos in merging clusters offer excel- 
lent opportunities for such comparisons. In particular, Markevitch et al. (2000) found two 
sharp X-ray surface brightness features in the merging cluster Abell 2142. They argue that 
these edges are formed by ram-pressure stripping during an ongoing merger in which the two 
components have passed each other once. Optical data indicate that the merger velocity is 
M ~ 1.5 and that cXc ~ 1200 km s"^ QOegerle, Hill, fc Fitchett 1995|) . The smaller com- 



ponent has i?gas ~ 50/i ^ kpc. Given the observed parameters of the cluster gas ([Henry fc 



Priel 1996|) , equation ([T7|) suggests that the velocity dispersion of the smaller component is 
ag ~ 700 km s^^. Our FDM model therefore predicts a mass enhancement (5S/S ~ 11% be- 
hind the smaller subcluster. In the CDM model, we predict a more gradual density increase 
with a peak amplitude 5S/S ~ 45%, although the system is at the limit of the validity of 
linear theory. Brightness edges similar to these have been observed in several other clusters 
( |Vikhlinin, Markevitch, fc Murray 2001| ; [Sun et al. 2001| ), indicating that such events are 



not uncommon. 

An alternate approach to tracing the mass distribution around the subhalo is to map the 
distribution of nearby dwarf galaxies (Markevitch & Vikhlinin 2001, private communication). 
Dwarfs behave as coUisionless test particles and so form a wake similar in structure to that 
shown in the top right panel of Figure 1. The amplitude of the variation in surface density of 
dwarfs near galaxies and groups can then be used to constrain the mass of the halo and hence 
the physics of the dark matter. Asymmetric distributions of dwarfs near infalling galaxies 
have been observed in the Coma cluster ([Conselice fc Gallagher 1998[ , 1999). However, the 
observed amplitudes of the density variations are 5S/S ~ 5 — 6, two orders of magnitude 
larger than our prediction for a galaxy moving through a rich cluster, so gravitational wakes 
are responsible for only a minor part of this enhancement. In fact, the observed excess of 
dwarfs in each system is contained near or within the tidal radius of the primary galaxy, 
indicating that the dwarfs may be bound to the primary or may have even been stripped from 
the primary due to ram pressure or strong tidal interactions with nearby galaxies ( [Conselice 



[fc Gallagher 199"^ ). Achieving observable signal-to-noise ratios on larger scales will most 
likely require targeted observations of the wakes of groups falling into clusters. 

For an impact parameter b outside the cluster core, S oc (1 + 6^/r^)~^'^, while 5S oc 
Ron{b) oc (1 -|- 6^/r^)~^/^. Therefore the relative signal is roughly independent of position 
throughout the cluster, although the absolute mass density (and hence the induced conver- 
gence and shear in lenses) will decrease as b increases. The increased number of cluster 
galaxies outside the core may make the envelopes more amenable to a statistical analysis. 
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4.2. Surface Brightness Maps 

Even without mass maps, it is possible to use observations of the ICM gas wake to 
constrain the dark matter properties. The ICM gas is best observed with X-rays, in which 
the emitted surface brightness is S* = /(e///47r)(i£. Here e// oc p^T^^"^ is the bremsstrahlung 
(free-free) volume emissivity. We assume in the following that the temperature of the gas 
varies adiabatically with its density (T oc p^^^), because recent studies have shown that 
heat conduction is suppressed in the ICM ( Vikhlinin, Markevitch, Forman, fc Jones 2001| ). 



Holding the gas isothermal decreases the signal by a small amount (see below) but does not 
significantly affect the structure of the observed feature. 

In searching for the wake, one may look for an abrupt drop in surface brightness at the 
edge of its cone. An estimation of the signal strength in this case is less straightforward 
analytically, because the surface brightness of the Mach cone produced by a point mass 
diverges (due to the singularity at the edge of the cone). Fortunately, the singularity is 
smoothed out by an extended perturber, and far from the perturber Agas <C 1 throughout 
the cone. In this region, we can expand the (1 + Agas)" that appears in the emissivity and 
retain only the first-order term. Here n = 7/3 if the gas temperature varies adiabatically 
with density and n = 2 if the gas is isothermal. We then find that 
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with the same assumptions as in equation (^). Here F{x, y, z, p) is the Gauss hypergeometric 
function (see |Gradshteyn fc Ryzhik 2000| , §9.10). The numerical factor is ~ 4 (or 5) if n = 2 



(or 7/3). As with the mass excess, the amplitude of the CDM wake is enhanced over the 
FDM wake by a factor ~ ^fdm- 

Equation (^) shows that the best targets for observations are the same as for lensing, 
either galaxies in groups or mergers of rich groups with clusters. If the dark matter cross- 
section increases with decreasing particle velocity as recently suggested ( |Hogan fc Dalcanton] 
20001 ; [Wyithe, 'lurner, fc Spergel 2000| ), galaxy groups may allow us to probe a more strongly- 



interacting regime that is closer to the FDM model. 

The bottom two panels in Figure 1 show the predicted surface brightness of the merging 
group and cluster described in the previous subsection. The wake from a halo moving through 
FDM (or CDM) is shown on the left (or right). The horizontal dashed lines again mark the 
points at which the central Mach cone intersects the cluster core. Note that the density 
enhancement is in fact greater than unity for < 2; < 1.5 and within the central Mach cone 
in the CDM case; however, the linear approximation remains valid in the remainder of this 
panel and in the left panel. 
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Figure 2 shows analogous surface brightness maps for galaxies moving through a group 
with (Jc = 500 km s^^, Uc = 2 x 10"^ cm~^, and Vc = 0.3 Mpc: the left panels show the wake 
for a galaxy similar to the Milky Way, with o"^ = 170 km s~^, and the right panels show the 
wake from a typical massive elliptical, with ag = 200 km s~^. Both galaxies have Ai = 1.5. 
Results for the FDM model are shown in the top panels while those for the CDM model 
are shown in the bottom panels. Horizontal dashed lines show where the central Mach cone 
intersects the edge of the group core. Note that the maps do not include emission from the 
hydrodynamic wake, stripped gas, or gas still bound to the perturber. We have also ignored 
any possible separation between the dark and gaseous halos (see §3.2). 

The sharp gradients near the central Mach cone in the CDM panels are not due to 
shocks (which we do not include in our formalism) but are instead due to the assumption 
of a singular core for the perturber. In a realistic case, the perturber will have a finite core, 
which smooths the gradient over the diameter of the core. To illustrate this effect, we show 
in Figure 3 the wake due to a uniform density perturber in the FDM (left panel) and CDM 
(right panel) models. The parameters are chosen identically to the right panels of Figure 
2, with the masses normalized to that of the corresponding isothermal sphere case. The 
smoothing of the central gradient by an extended core should help distinguish CDM wakes 
from the sharp cones expected in the FDM model. The morphological differences between 
the singular and smoothed cores are much less dramatic in the FDM case because the halos 
are truncated at much smaller radii in that case. 

Limits on the nature of the dark matter from surface brightness maps are less direct 
than those from lensing maps, because we must use the ICM gas wake to infer the size of 
the dark matter halo. The first observational constraint comes from the width of the edge of 
the gravitational Mach cone, which is in turn determined by Rg. Let the observable size of 
the gaseous halo be -Rgas- For FDM, we expect that Rg ~ i?gas so that the edge of the cone 
is tight and follows closely the observable edge of the perturber. If, on the other hand, the 
dark matter is collisionless, Rg ~ a~^gRgg,s 3> -Rgas, and the edge of the gravitational cone is 
very extended. However, because the ICM gas is collisional, we expect hydrodynamic and 
stripping wakes to be present as well. The characteristic widths of these wakes is ~ -Rgas; so 
we expect a well-defined cone to trail behind the gaseous halo even in the CDM model. 

Fortunately, there remain signatures in the surface brightness distribution that can dif- 
ferentiate between the FDM and CDM models. First, stripped gas can be observationally 
distinguished from ICM gas through measurements of metallicity or gas temperature (the 
stripped gas will presumably be cooler and more metal- rich than the ICM gas). The hydro- 
dynamic wake, like the gravitational wake, is composed of perturbed ICM gas, which will, if 
anything, be heated by compression in the wake. Therefore we can separate these two wakes 
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Fig. 2. — X-ray surface brightness maps for galaxies moving with M. = 1.5 through a 
background cluster with cTc = 500 km s~^, ric 



2 X 10 ^ cm ^ and r, 



(Tg = 170 km s ^, FDM model. To|? ng'/it.- o"„ = 200 km s 



0.3 Mpc. Top left: 
FDM model. Bottom left: 



(jg = 170 km s 



CDM model. Bottom right: ag = 200 km s~^, CDM model. Contours 
are set at 5% of the peak surface brightness of the background cluster in each panel. The 
merger occurs in the plane of the sky with the perturber moving along the z-axis. Axes 
(in cylindrical coordinates) are in units of Rq appropriate to each panel (see Eq. |T9[] ). 
Horizontal dashed lines mark the intersection of the Mach cone with the cluster core. 
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Fig. 3. — X-ray surface brightness maps for galaxies moving with M. = 1.5 through a 
background cluster with cTc = 500 km s~^, ric = 2 x 10~^ cm~^ and Vc = 0.3 Mpc. The 
perturber is assumed to have a uniform density, with total mass normalized to a singular 
isothermal sphere of CTg = 200 km s~^. Left: FDM model. Right: CDM model. Contours 
are set at 5% of the peak surface brightness of the background cluster. The merger occurs 
in the plane of the sky with the perturber moving along the 2;-axis. Axes (in cylindrical 
coordinates) are in units of Rq appropriate to each panel (see Eq. |]19|). Horizontal dashed 
lines mark the intersection of the Mach cone with the cluster core. 
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from contamination by the stripped gas. 

Even without such a separation, a robust indicator of the CDM model is a surface 
brightness enhancement leading the hydrodjTiamic bow shock. Because FDM halos are the 
same size as the gas halo, causality dictates that their effects be confined behind the bow 
shock leading this body. The distance between the bow shock and the front surface of the 
body may be estimated if the Mach number and geometry are known (|Vikhlinin, Markevitch 



|fc Murray 200 1| ; [Schreier 1982|) . Any enhancement leading this shock must be due to the 



gravitational effect of a CDM halo. Furthermore, in the CDM model, the gas halo may trail 
the dark matter halo because the latter does not experience ram pressure. In this case we 
would observe the gravitational brightness enhancement to be displaced even farther in front 
of the hydrodynamic Mach cone. 

Abell 2142 is again a prime candidate for observations. Our model predicts a sur- 
face brightness excess 6S/S ~ 0.4-0.5 in an FDM cone behind the object and a peak 
excess 6S/S > 1 in the CDM model, which should be easily observable in deep observa- 
tions with Chandra. Galactic candidates also exist. NGC 1404, an elliptical galaxy with 
ag ~ 150 km s~^ ( |Winsall fc Freeman 1993| ) appears to be falling into the core of the For- 



nax group, which has 0"c ~ 370 km s ^ ( Drinkwater, Gregg, fc CoUess 200 1| ). The infall 



velocity of the system is 7W ^ 1.1 ( [Ferguson 1989|) . Using the X-ray data of Jones et al. 



(1997), we predict SS/S ~ 0.12-0.15(rc/125 kpc) in the FDM model, where r^ is the assumed 
core radius of the group. The peak brightness excess in the CDM model is expected to be 
6S/S ~ 0.6(rc/125 kpc). Paolillo et al. (2001) have made deep ROSAT observations of this 
system, and they see evidence for a small plume behind the galaxy. They do not see evidence 
for any extended cone emission, although the level we predict is equal to only one of their 
contours. This system therefore provides some weak evidence against the FDM model. The 
morphology is more consistent with a CDM wake, especially because there appears to be 
substantial emission leading the core of the galaxy. However, higher resolution observations 
are necessary in order to precisely determine the ram pressure edge of the galaxy. Further- 
more, the trailing plume is not as extended as one would expect for the CDM case, although 
projection effects may disguise the wake. 



4.3. Comparison to Numerical Simulations 



Numerous groups have simulated cluster mergers (e.g., [Roettiger, Loken, fc Burns 1997 



[Ricker 1998| ; [Ricker fc Sarazin 2001| ; [Ritchie fc Thomas 2001|) and galaxy/ICM interactions 



e.g., [Stevens et al. 199^ ; [Quilis et al. 2000| ) over the past several years. Unfortunately, 



neither of these classes of simulations is ideal for analysis of the gravitational wakes. In the 
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latter case, the simulation box size is on the order of the galactic diameter, where nonlinear 
effects are strongest and where stripped material dominates the wake. In the former case, 
the mergers studied are generally in the nonlinear regime (most commonly, two equal mass 
clusters collide) and in complicated environments from which it is difficult to isolate the 
effects of the gravitational wake. 

In fact, our formalism does not apply to existing studies of merging clusters in the CDM 
model because the background density varies rapidly. Such simulations ([Roettiger, Loken, fc 



Burns 1997| ; [Ricker fc Sarazin 2001| ; [Ritchie fc Thomas 200 1| ) have re ~ 1, so our formalism 



fails in all of these cases, even if linear theory would otherwise be valid (see equation [0). In 
order to test our model and evaluate the importance of shocks and stripped gas, we suggest 
simulations with relatively uniform backgrounds fulfilling the criterion fc ^ 1. 

Recent simulations also show that shocks become significantly less important as the 
mass of the perturber decreases because the increase in the entropy of the core gas is small 
for such mergers ( [Ritchie fc 'I'homas 2001| ). Roettiger, Loken, & Burns (1997) also find that 
gradients in the Mach number of the gas decrease as the mass ratio of the merger decreases, 
indicating again that shocks become significantly weaker. We therefore expect linear theory 
to be an adequate approximation in the limit of a small perturber mass. 

To date, there has been little investigation of the direct effects of cluster mergers in 
the FDM regime. Ricker (1998) simulated the merger of two equally massive gaseous halos, 
but we do not expect linear theory to apply in this case. We therefore cannot use existing 
simulations to test the accuracy of our predictions for the FDM model. 

A complementary set of simulations to those discussed above have been performed by 
Sanchez-Salcedo & Brandenburg (1999, 2001), who investigated the accuracy of the dynam- 
ical friction force found by Ostriker (1999) for a point source traveling through a collisional 
medium (equation ^). They do not directly compare the predicted and simulated wake 
structures, but they do find that the drag force expected from linear theory is accurate to 
within several percent of the simulation results, provided that the Coulomb logarithm is cho- 
sen appropriately ( [Sanchez-Salcedo fc Brandenburg 1999| ). They have also studied the drag 
force on a perturber in a spherically symmetric Plummer model potential ( [S anchez- S alcedo 



fe Brandenburg 2001[). In their model, Vc ~ 6, so we would expect our formalism to be 



marginally valid. Again, they do not analyze the wake structure, but they do find that the 
drag force on the perturber is approximately that given by linear theory if the characteristics 
of the background medium at the instantaneous location of the perturber are used. However, 
they do not include the hydrodynamic effects discussed in §2.3, so we cannot evaluate their 
importance relative to gravitational effects. 
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5. Discussion 

We have calculated the gravitational density wake of a perturber moving through a 
uniform collisional or coUisionless medium. We also described how the dark matter of a 
galaxy or group moving through a cluster will be stripped in the cases of fluid-like dark 
matter (FDM) and coUisionless dark matter (CDM): the former experiences ram pressure 
stripping, while the latter undergoes tidal stripping only. Finally, we calculated the effects 
of the wakes in the two different regimes on both the projected cluster mass and the X-ray 
surface brightness. 

In the FDM model, a supersonic perturber generates a dark matter wake with a well- 
defined Mach cone in which the surface density increases substantially in a relatively narrow 
region; in the CDM model, the surface density increases more slowly over a broad region but 
is of larger amplitude. Because the intracluster gas is collisional, we expect the gravitational 
wake in the gas to form a Mach cone in both the FDM and CDM models. Nevertheless, 
because FDM halos are truncated at much smaller radii, the wakes in the two models have 
significant morphological differences that can be used to distinguish them. 

We found that the signals due to these wakes are potentially observable if the perturber 
moves supersonically. For a large elliptical galaxy (o"g ~ 250 km s~^) in a rich cluster core 
(cTc ~ 1000 km s~^, ?2c ~ 5 X 10~^ cm~^, Tc ~ 0.4 Mpc), the variation in the surface density 
is (5S/S ~ 0.7% across the edge of the FDM Mach cone, while the variation in surface 
brightness is 6S/S ~ 3-4% for the FDM model. In the CDM model, we expect a more 
diffuse wake with signal strength ~ 6 times those in the FDM model. Fortunately, both 
signals increase rapidly with the relative velocity dispersions of the perturber and cluster 
(see equation ||2^). Therefore, the most favorable systems for observations have a large 
Ogjoc- either a gas-rich group moving through a cluster or a galaxy moving in a group. 
For example, the expected signals for a typical galaxy [pg ~ 170 km s~^) in a galaxy group 
(o-c ~ 500 km s-\ nc ~ 2 X 10"=^ cm'^, r, ~ 0.3 Mpc) are (5S/E ~ 1.3% and bSjS ~ 7% in 
the FDM model (or 5S/S ~ 6% and bSjS ~ 60% in the CDM model). Statistical methods 
may be used to enhance the signal-to-noise ratio; a related technique, used on ROSAT data 
to search for small-scale wakes due to ram pressure stripping, has been developed by Drake 
et al. (2000). The relative signal strength is insensitive to the location within the cluster, so 
it may be possible to perform such a statistical analysis over a large volume with sufficiently 
deep imaging. 

The resolution and sensitivity of the Chandra satellite allows detailed study of the inter- 
actions between individual subhalos and the ICM. First, surface brightness and temperature 
observations on both sides of the bow shock leading the perturber (caused by hydrodynamic 
effects unrelated to the gravitational wake; see §2.3) yield the Mach number of the perturber 
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through the shock jump conditions (e.g., [Landau fc Litshitz 1959 ). We can combine this 



information with the observed radial velocity of the perturber or with the opening angle of 
the shock at large distances to find the inclination angle of the cone. When combined with 
optical observations (yielding ag and Uc) and X-ray observations (yielding ric, Tc, and gas 
temperature), we can make an unambiguous prediction of the enhancement expected from 
the gravitational wake for the FDM and CDM models. This does not include the enhance- 
ment due to stripped gas or deflection of the ICM particles, but for each model it gives a 
lower limit to the expected enhancement. Therefore, observations showing enhancements 
smaller than those predicted by the CDM model can be used to falsify that model. In addi- 
tion, temperature maps differentiating the cool stripped gas from hot ICM gas can help to 
constrain the cause of any particular observed enhancement. 

Existing X-ray maps of cluster mergers and galaxy motion in groups have shown no 
evidence for the sharply defined Mach cones that would be expected in the FDM model 
(e.g., Paolillo et al. 200l| ); however, there is very little existing data of a high enough quality 



to address this question. The new generation of X-ray satellites have the resolution and 
sensitivity to observe these features with deep observations. 

Our model predicts two systematic variations in the wakes. First, the extent of the 
galaxy halos should increase with distance from the cluster center: both tidal truncation 
(for CDM halos) and ram pressure stripping (for FDM halos) become more effective as the 
ambient density increases. We expect Rg oc n~^/^ in both cases, so this systematic variation 
cannot be used to test the nature of dark matter, but it can be used as a consistency check 
on statistical results for halo sizes. 

The second systematic trend can be used to probe the coUisional nature of the dark 
matter. Tidal truncation depends only on the relative gravitational potentials of the per- 
turber and background cluster and is thus independent of the velocity of the perturber. Ram 
pressure, on the other hand, increases with velocity. We therefore expect that halo size will 
decrease with increasing peculiar velocity relative to the cluster in the FDM case (where 
Rg oc TW^^; see equations Jl^ and ||18[) but not in the CDM case. 

Our model for FDM assumes that the dark matter is a perfect fluid, while many SIDM 
models treat the dark matter as having a mean-free-path of order the size of the system. 
In such cases, the structure of the wake (and hence its observability) will lay somewhere 
between the extreme cases described here. Unfortunately, these more moderate cases are 
not amenable to analytic treatment. Instead, one must make use of Monte Carlo techniques 
in numerical simulations (e.g., [Burkert 20001 ; [Kochanek fc White 2000| ) to analyze the wake 



structure. Such simulations would also help to address some of the other shortcomings of 
our model (e.g., the use of linear perturbation theory in the collisional case, the assumption 
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of a uniform density background, the structure of the hydrodynamic wake, and the role of 
shocks; see also §4.3). 

Finally, we outlined in §3 a number of consequences of the different truncation mecha- 
nisms acting on halos in clusters, including stellar mass loss and its effects on rotation curves, 
the fundamental plane, and local group dwarfs, relaxation times of galaxies in clusters, and 
separation of the dark and gaseous components of galaxy groups. While indirect, these ef- 
fects do not rely on the fluid approximation and may therefore be useful in constraining the 
dark matter properties between the fluid and collisionless limits. 

We thank W. Forman, C Jones, M. Markevitch, and A. Vikhlinin for very helpful 
comments on the manuscript. We would also like to thank P. Ricker for making unpublished 
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5-7039, 5-7768, and NSF grants AST-9900877, AST-0071019 for AL. SRF acknowledges the 
support of an NSF graduate fellowship. 
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